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A family of wave solutions to the full Einstein equations in AdSs geometry is derived. These 
background solutions give by duality the response of M— 4 SYM at strong coupling to an arbitrary 
distribution of fast moving external sources for the energy-momentum tensor operator. We discuss 
the similarities of these solutions with the Color Glass Condensate effective theory for QCD at high 
energy. 
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I. INTRODUCTION 

The AdS/CFT correspondence [1-3] and other 
gauge/gravity dualities have brought new insights into 
strongly coupled gauge theories. Since the data have 
provided evidences that the hot and dense matter pro- 
duced in heavy ion collisions at RHIC should be strongly 
coupled [4], AdS/CFT has become an interesting tool 
in order to explore subjects related to RHIC physics. 
The quark-gluon plasma (QGP) presumably produced 
at RHIC is deconfined and at least approximately locally 
thermalized. In this regime, QCD acquires some similar- 
ities with Af— 4 SYM at finite temperature, which can 
justify the use of AdS/CFT. In this context, the hydrody- 
namical behavior of the long wavelength modes of Af= 4 
SYM at finite temperature has been understood and re- 
lated to black brane physics [5, 6]. This has allowed to 
calculate the shear viscosity [7] and other transport co- 
efficients [8-10] of a 7V= 4 SYM plasma. Assuming the 
symmetries of the Bjorkcn flow, it has been shown that 
hydrodynamics is the only consistent behavior at late 
proper time for the strongly coupled plasma [11-13]. The 
picture of the plasma emerging from these studies is in 
qualitative agreement with the one extracted from the 
data. 

Gkion saturation [14] dominates the dynamics of the 
nucleus before the collision, and determines the initial 
conditions for the medium produced. The correct frame- 
work allowing to address that phenomenon from QCD 
first principles is the Color Glass Condensate effective 
theory (CGC) (see [15, 16] for reviews), or equivalent 
formalisms, such as the one developed by Balitsky [17]. 
The coupling is weak at high enough energy, due to the 
properties of gluon saturation. Among other successes, 
the CGC provided a correct prediction for d-Au collisions 
at RHIC [18], where initial condition effects are essential. 

One of the main theoretical problems about heavy 
ion collisions is to understand the thermalization pro- 
cess leading from the weakly coupled out-of-equilibrium 
Glasma [19], produced by the colhsion of two sheets of 
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CGC, to the strongly coupled and thermalized QGP. Let 
us consider the possibility that thermalization is mainly 
a strong coupling phenomenon. Then, thermalization 
should be addressed by building the gravity dual of an 
expanding medium, starting from very short proper time. 
In that regime, the physics is determined by the proper- 
ties of the nucleus before the collision. Hence, it would 
be useful to construct the gravity dual of a fast moving 
object mimicking the physics of the CGC. 



The obvious problem is that, as a conformal theory, 
J\f= A SYM admits no states analogous to nucleus or 
hadrons. There are two main strategies allowing to avoid 
that difficulty. The first strategy, which is commonly 
used, is to modify the gravity background dual to J\f= 4 
SYM, in order to obtain confinement in the IR. It is how- 
ever difficult to maintain a gauge / gravity correspondence 
in that context. Less drastic modifications of the grav- 
ity setup such as the introduction of point-like sources 
in the bulk [20, 21] seem to provide models for a fast 
moving nucleus in the boundary theory. However, the 
modification of the boundary field theory corresponding 
to the addition of sources in the bulk is not well under- 
stood. The second possible strategy, developed in the 
present paper, is the following. We restrict ourselves to 
the original AdS/CFT correspondence [1-3] , keeping con- 
sistently non-zero sources of J\f = 4 SYM, and we make 
use of the results of holographic renormalization [22, 23]. 
Suitably chosen distributions of sources generate Af = 4 
SYM fields whose energy-momentum tensor is a reason- 
able approximation to the one of a fast moving nucleus in 
QCD. Such setup is well defined both on the field theory 
side and on the gravity side of the correspondence, and 
has some similarities with the weakly coupled CGC. 



In section II, we derive a family of background solu- 
tions of Einstein equations. In section HI, the boundary 
field theory counterparts of these solutions are given, re- 
lying on the results of holographic renormalization. In 
section IV, we discuss the similarities and differences of 
our solutions with the CGC formalism, and suggest how 
to model a highly boosted nucleus at strong coupling. 
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II. A FAMILY OF EXACT SOLUTIONS OF 
EINSTEIN EQUATIONS 

Let us look for gravity duals of Af~ 4 SYM with fast 
moving matter in the (+) direction on the light-cone. 
The corresponding metric writes as the AdS^ metric plus 
an additional part to be found. In the high energy limit, 
we keep for that new part only the component which 
is enhanced twice by Lorentz factors with respect to the 
rest frame of the matter distribution. Hence, generalizing 
the Shockwave solution proposed in [11], we consider the 
Ansatz^ for the metric Gij, in light-cone coordinates. 



dz^ —2dx'^ dx -|-dxj_^-|-6(a;"'", a; ,xj_,2;)di 



(1) 

where xj^ is a vector of the 2-dimensional transverse 
plane. We do not modify other supergravity background 
fields, so that the cosmological constant should remain 
A = — Hence, one can write the Einstein equations 
as 



(2) 



The {z, — ), {i, — ) and (+, — ) components of the Einstein 
equations for our Ansatz (1) reduce respectively to 

d,d+b = 0, d,d+b = and dlb = . (3) 

In the function b{x^ z), the linear term in x^ 
allowed by the equations (3) does not encode interesting 
physics of the boundary field theory. Hence, we will re- 
strict ourselves, in the following, to functions b{x~ ,x±,z) 
independent of a;+. 

Now, the non-vanishing Christoffel symbols for our 
Ansatz are 

ri_ = dzb{x~ ,yi±,z) + - b{x~,xj_,z) 
I z 



1 



d^b{x ,x_L,2;) 
1 



rl^ = r+_ = dib{x-,x^,z).. 



(4) 



in addition to the Christoffel symbols of pure AdS^ ge- 
ometry, which turn out to stay unmodified^ 



n_=ri+ = rL = r^ 



zl 



(5) 



^ We write with capital letters I,J,... 5-dimensionaI AdSs-likc in- 
dices, with Greek letters fi,u,... the 4-dimensional Minkowski- 
like indices, and with small Latin letters i,j,... indices for the 
2-diuKniKional transverse plane. 

^ There is no summation over repeated indices in formula (5). 



Calculating the Ricci curvature TZij, one finds that all 
the Einstein equations are trivially satisfied, except the 
one corresponding to the (— , — ) component, which gives 

3 

d'^b{x~ ,x±,z) dzb{x^ ,xj_, z) + Aj_b{x~ ,xj_, z) = 0. 

(6) 

In many studies of gauge/gravity duality, one considers 
small perturbations on the top of a known background, 
in order to linearize the problem. It should be remarked 
that here, by contrast, the perturbation b can be large, 
as the exact Einstein equations turn out to be linear in b, 
thanks to the peculiar structure of the Ansatz (1). Per- 
forming a Fourier transform with respect to the trans- 
verse position x_L , one gets 

3 ~ 

d^b{x~ jk.^, z) dzb{x~ yli^, z) — b{x~ ,)s.±, z) = . 

z 

(7) 

For k^ ^ 0, the generic solutions of equation (7) write 



b{x ,kj_,2) = ci{x ,kj_) 



2 ^2 



kiz 



K2{k^z) 



+ C2{x- kiz' hik^z) , (8) 

where I2 and K2 arc the modified Bcsscl functions of first 
and second kind, k±_ = |k_L|, and Ci and 62 are two arbi- 
trary functions. z''K2{k_\_z) decays exponentially to zero 
as z — >■ 00 whereas 2:^/2 (fc_Lz) blows up exponentially. 
Our Poincare coordinates system covers only half of the 
AdS^ space. For that reason the points with z ^ 00 are 
in the interior of the manifold. Hence, only z'^K2{k^z) 
is regular in the bulk. Since the function b disappears 
from the calculation of the square of the Riemann tensor 
and any other scalar quantities built from the Riemann 
curvature, one could think that our Ansatz never leads 
to curvature singularities. However, a similar family of 
solutions - up to dimension and signature - has been an- 
alyzed [24], and the singularity 2 00, if present, has 
been shown to be a p. p. curvature singularity, as defined 
in [25]. This remark should hold in our case, meaning 
that the z''l2{k±z) solutions lead to a naked curvature 
singularity at z — > 00. 

z'l2{k±z) might be acceptable only as the small z be- 
havior of a solution with a particular additional object in 
the bulk, mapping to a modification of the IR behavior 
of the boundary field theory. In the present study, we 
will not discuss further that possibility, and keep only 
the regular modes with z''K2{kxz). 

The solutions of equation (7) for the k_L = modes 
writes 



b{x ,\s.± = 0, z) = di{x ) + z'^ d2{x ). 



(9) 



As for the generic case, we will discard the second term, 
which should lead to a naked curvature singularity at 
z — >■ 00. One should notice that this singular mode is 
precisely the Shockwave solution proposed in [11] as a 
model for a highly boosted nucleus, and used in [26, 27] in 
order to model nucleus- nucleus collisions or deep inelastic 
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scattering on a nucleus. However, one should realize that 
such a solution which is singular in the bulk does not map 
to a solution of pure N=A SYM in the vacuum. 



As 



K2{k_ 



1 for z ^ 0, we merge the two 



cases, and write our solution for arbitrary kj_ as 



h{x ,k_L,z) = ci{x ,k_L) 



U2 ,2 



(10) 



with c\{x , k j_ = 0) = di (a; ). We have a family of reg- 
ular background solutions (1,10) of Einstein equations. 

In that family, the only solutions which are homo- 
geneous in the xx transverse plane are of the type 

6(x~,x_L,2:) = h{x~). In that case, performing the 

change of variable x+ i"*" = .t+ — dub{u)/2, one 
obtains the undeformed AdS^ space, with i+ playing the 
role of Hence, these homogeneous solutions seems 
not to encode any no n- trivial information about J\f= 4 
SYM. We will come back on that issue later on. 

Lot us now discuss the content of the solutions (1,10) 
in terms of the dual field theory. 



(13) is3 



(T,^) = [Tr(5(2)ff(o;5(2)) - (TV 5(2))^ 



9{0h 



+5(4);^^-^ (5(2)5(0)5(2) ),^i. + ^(Tr g{2))g{2)^iu \ (14) 



up to renormalization-scheme dependent terms propor- 
tional to /l(4)^i^. 

If one imposes to have the Minkowski metric on the 
boundary 5(o)/Lti/ = Tlnv once the functional derivation in 
equation (13) is done, one calculates (T^j/) in the vacuum. 
Hence, on the field theory side, 5(o)^;y — Vt^v corresponds 
to a distribution of sources of T^i/. The gravitational 
background solutions (1,10) are thus dual to the JV= 4 
SYM theory with non-vanishing external sources J^,^ of 
the energy-momentum tensor, given by 

J^^(a;",xj_) = ,9(o),„Xa;-,x_L,2: = 0) - r?^^ 

= dn-6„- b{x~,x±,z = 0). (15) 



III. DUAL FIELD THEORETICAL CONTENT 
FROM HOLOGRAPHIC RENORMALIZATION 

The asymptotic behavior near the boundary of back- 
ground metrics of the type 



dz^ + g^i^{x'',z)dx''dx'' 



(11) 



has been analyzed in [22]. When the z = Q boundary is 
4-dimensional, one finds the expansion 

g^u{x^, z) = g(Q)^,v{x^) + z^ g{i)tLv{x'') z^ g(4)^iu{x'') 
+z^logiz^) /i(4)^.(a;^) + O (zHogz) ,(12) 



That corresponds to a distribution of sources moving to 
the right at the speed of light, with arbitrary longitudinal 
and transverse profiles. 

In the case of the solutions (1,10), the only non- 
vanishing component of 5(2)^1^, 5(4)/^;/ and /i(4)^y is the 
(-,-) one. Moreover one has g'^^J = 0. Therefore, all 
traces and other nonlinear terms in the metric coefficients 
present in (14) and in other formulae given in [22] vanish, 
which leads to 



L3 



47rG 



JV 



5(4)- 



^4)- 



5(2)- 



8 ^^5(2)- 
A_l5(0)— 



(16) 
(17) 
(18) 



where 5(2)/ii/ and /i(4)^i/ are completely determined by 
5(0) MJ^- "^^^ expectation value of the energy- momentum 
tensor (T^,y) of the dual conformal field theory living 
on the boundary is calculated according to the usual 
AdS/CFT prescription [2, 3] and with the appropriate 
holographic renormalization [22, 23] as 



^-det 3(0) Sg 



■JW {Sgr,reg[g{0)\ - Sgr,ct[g(Q)\) , 



(0) 



(13) 

where Sgr.reg is a regularized version of the on-shell grav- 
itational action, and Sgr.ct contains the necessary coun- 
terterms. Regularization and renormalization are re- 
quired due to IR divergences near the boundary in the 
gravity theory which map to UV divergences in the dual 
field theory. The general result [22] of the calculation 



Thanks to the expansion* 



P 72 

K^ik^z) = 1 



klz^ 



3 „ 
-2-27. 



~32" 



-2 log 



k±z 



+ O{z''logz), (19) 



one can identify the transverse Fourier transform of the 
coefficients of the expansion of the metric (1,10) near the 



^ The traces in formula (14) are taken with respect to the metric 
^ — 0.577 is the Euler-Mascheroni constant. 
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boundary, as 

5(0)— (a;",k_L) = ci(a;",k_L) 

5(2)— (a;",k_L) 

.(x-,kx) + log(22) /,, 



^2 

ci(a;-,kx) 



ff(4)- 
32 



-21og('%^ 



(4)- 



(20) 
(21) 

ci(a;-,ki). (22) 



From the expressions (20), (21) and (22), one checks that 
the relations (17) and (18) are satisfied, and thus the cor- 
responding generic relations found in [22] are also satis- 
fied, as they should. As z is dimensionful, the left hand 
side of (22) is a formal and ambiguous expression. One 
has to introduce an arbitrary momentum scale in or- 
der to write log {k±z/2) = log {k±/fiR) + log in 
the right hand side of (22). The ambiguity in the choice 
of reflects the ambiguity in the counterterms used in 
the formula (13), and thus the choice of the renormal- 
ization scale in the boundary field theory. Moreover, it 
is possible to split arbitrarily the constant terms in the 

brackets in equation (22) into a contribution to 5(4) 

and a contribution to log(z^) /i(4) . This last prop- 
erty reflects the freedom in the choice of renormalization 
scheme in the boundary field theory. Choosing for sim- 
plicity to put as many terms as possible in log(^;^) /i(4) , 

one finds the Fourier transform of the only non vanishing 
component of (T^,^) 



using the AdS/CFT dictionary to express the 5- 
dimensional Newton constant as Gat = One should 

remark that (T^) = 0. Hence, the presence of the fast 
moving external sources does not bring conformal anoma- 
lies for A/'=4 SYM. 

Due to the explicit presence of sources on the gauge 
theory side, our construction can be understood as an 
operator deformation of the AdS/CFT correspondence, 
by the T operator of J\f= 4 SYM. Since the metric per- 
turbation b{x~ , xj^, 2;) has a finite limit on the boundary, 
this deformation is an exactly marginal one. Due to its 
dimension 4, T^,y is indeed a marginal operator of J\f— 4 
SYM. Hence, in the ultraviolet on the gauge theory side, 
we do not have the pure J\f= 4 SYM theory but exactly 
marginal deformations of it. 

Before closing that section, let us come back to the 
homogeneous solutions b{x~ ,x±,z) = b{x~). From 
the general results of [22], one finds easily that these 
solutions correspond in the dual gauge theory to a 

distribution of sources J (x~) which do not induce any 

non-trivial (T )(a;~), confirming the triviality of these 

homogeneous solutions. Intuitively, one can understand 
that from the M= 4 SYM side as follows. After a very 



strong boost along x'^ of a generic field configuration, the 
dominant components of the N= A SYM field strength 
should bc! the ones, and the dominant component 
of the (r^„) should bc (T__) oc (-F^,)^ In the case of a 

homogeneous distribution of sources J (a^ ), there is 

no preferred direction in the transverse plane, and thus 
has to vanish, and (T--) as well. By contrast, if 

J (a;~,x^) depend on xj_, the transverse gradient of 

J gives a preferred direction in the transverse plane, 

allowing a non- vanishing (T ). 



IV. DISCUSSION: STRONGLY COUPLED 
ANALOG OF A COLOR GLASS CONDENSATE 

Our background solutions (1,10) are dual to A/'=4 SYM 
with rightmoving sources on the lightcone. The formula 
(23) gives the energy and momentum of the SYM fields 
induced by the sources, including the potential energy 
of the fields due to the presence of the sources. This 
picture is reminiscent of the CGC framework, where a 
highly boosted nucleus is described by fast moving color 
sources representing the long-lived partons and by the 
quasi-classical gluon field F"_ they induce. The only 
non-zero component of (Tjuj/) for this gluon field is also 
(T ). An important difference between the two frame- 
works is that the CGC is naturally formulated in terms of 
colored fields whereas AdS/CFT gives access only to col- 
orless objects, which complicates a detailed comparison 
between the two setups. Supplemented by a distribu- 
tion of sources both peaked in the longitudinal direction 
and localized in a finite region of the transverse plane, 
the metric (1,10) give a dual model for a fast moving 
nucleus, which can be easily used to study high energy 
scattering at strong coupling, including effects of confine- 
ment on the geometry of the collision. Most observables 
usually considered are not sensitive to other, truly dy- 
namical, effects of confinement, which happens only d\ir- 
ing hadronization, and which are not implemented in the 
standard AdS/CFT approach. In principle, one can push 
further the analogy with the CGC framework by doing 
the calculation of an observable for a given J , and av- 
erage the result over J as done in the CGC, in order 

to take event-by-event fluctuations of the nucleus wave 
function into account. 

In summary, we have derived a large family of wave 
solutions (1,10) of the full Einstein equations in AdS^ 
geometry. Relying on the AdS/CFT correspondence 
and its holographic renormalization, we have shown that 
these gravity backgrounds are dual to M= 4 SYM with 
arbitrary distributions of rightmoving sources on the 
lightcone. This framework is a strongly coupled analog 
of the CGC, and provides a laboratory, suitable for the 
investigation of high energy scattering in A/'= 4 SYM 
and thermalization at strong coupling. 
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Note added: The paper [28] appeared on the arXiv 
a few months after the present manuscript. In [28], the 
authors discuss deep inelastic scattering on targets given 
holographically by various Shockwave backgrounds, in- 
cluding the ones derived in the present study. The au- 
thors of [28] discard the K2 solutions we are keeping (see 
Eq. (8)), and consider the I2 solutions we are discarding. 
However, there is not necessary a contradiction between 
the two choices, but possibly just two different prescrip- 
tions related to two different types of holographic heavy 
ion collisions. 

Indeed, the purpose of [28] is to find, using deep inelas- 
tic scattering, what kind of Shockwave gives the content 
expected for a nucleus in a confined version of TV = 4 
SYM at strong coupling. It is shown in [28] that the I2 
solutions are suitable in that case, provided one modifies 
the background at large z in order to implement confine- 
ment. From the studies of high energy scattering within 
AdS/CFT, one expects the collision of two such strongly 
coupled TV = 4 SYM nuclei to be quite different from 
real-world heavy ion collisions e.g. at RHIC, which show 
rather weak stopping properties. 

By contrast, the collision of two K2 Shockwaves should, 
at very early proper time, mimic more closely weakly cou- 
pled QCD heavy ion collisions, with possibly a Glasma 



[1] J. M. Maldacena, Adv. Theor. Math. Phys. 2, 231 (1998) 

[Int. J. Theor. Phys. 38, 1113 (1999)]. 
[2] S. S. Gubser, I. R. Klebanov and A. M. Polyakov, Phys. 

Lett. B 428, 105 (1998). 
[3] E. Wittcn, Adv. Theor. Math. Phys. 2, 253 (1998). 
[4] E. V. Shm-yak, Nucl. Phys. A 750, 64 (2005). 
[5] G. Policastro, D. T. Son and A. O. Starinets, 

JHEP 0209, 043 (2002), JHEP 0212, 054 (2002). 
[6] S. Bhattacharyya, V. E. Hubony, S. Minwalla and 

M. Rangamani, JHEP 0802, 045 (2008). 
[7] G. Policastro, D. T. Son and A. O. Starinets, Phys. Rev. 

Lett. 87, 081601 (2001). 
[8] M. P. Heller and R. A. Janik, Phys. Rev. D 76, 025027 

(2007) . 

[9] R. Baler, P. Romatschke, D. T. Son, A. O. Starinets and 
M. A. Stcphanov, JHEP 0804, 100 (2008). 
[10] M. Natsuume and T. Okamura, Phys. Rev. D 77, 066014 

(2008) [Erratum-ibid. D 78, 089902 (2008)]. 

[11] R. A. Janik and R. B. Peschanski, Phys. Rev. D 73, 
045013 (2006). 

[12] M. P. Holler, P. Surowka, R. Loganayagam, M. Spalinski 
and S. E. Vazquez, Phys. Rev. Lett. 102, 041601 (2009). 

[13] S. Kinoshita, S. Mukohyama, S. Nakamura and 
K. y Oda, Prog. Theor. Phys. 121, 121 (2009). 

[14] L. V. Gribov, E. M. Levin and M. G. Ryskin, Phys. Rept. 
100, 1 (1983). 



stage. Indeed, the presence of boundary sources allows 
to impose a high degree of transparency, and an approx- 
imate boost-invariance in the central region. However, 
these K2 Shockwaves presumably cannot be obtained 
from the ultrarelativistic limit of nuclear states of a con- 
fined strongly coupled gauge theory. 

Moreover, the problem of positive energy condition vi- 
olation found in [26, 27] for the collision of two z'^ shock- 
waves remains in the case of I2 Shockwaves, but should 
disappear in the case of K2 Shockwaves, according to the 
structure of Einstein equations. 

In summary, the K2 solutions maybe more suitable to 
build an holographic model for real-world QCD heavy 
ion collisions, although the I2 solutions are more closely 
related to boosted nuclei in a strongly coupled confined 
gauge theory. 



Acknowledgments 

I acknowledge losif Bena, Edmond lancu, Romuald 
Janik, Larry McLerran and Robi Peschanski for fruitful 
discussions and useful comments. 



[15] E. lancu, A. Leonidov and L. McLerran, arXivihep- 
ph/0202270. 

[16] E. lancu and R. Venugopalan, arXiv:hep-ph/0303204. 
[17] I. Balitsky arXiv:hep-ph/0101042. 

[18] D. Kharzccv, Y. V. Kovchcgov and K. Tuchin, Phys. Rev. 

D 68, 094013 (2003). 
[19] T. Lappi and L. McLerran, Nucl. Phys. A 772, 200 

(2006). 

[20] H. Nastasc, arXiv:hcp-th/0410124. 

[21] S. S. Gubser, S. S. Pufu and A. Yarom, Phys. Rev. D 78, 
066014 (2008). 

[22] S. de Haro, S. N. Solodukhin and K. Skenderis, Commun. 
Math. Phys. 217, 595 (2001). 

[23] M. Bianchi, D. Z. Freedman and K. Skenderis, Nucl. 
Phys. B 631, 159 (2002). 

[24] J. Podolsky Class. Quant. Grav. 15, 719 (1998). 

[25] S. W. Hawking and J. F. R. Ellis, The Large Scale Struc- 
ture of Space- Time (Cambridge: Cambridge University 
Press, 1973). 

[26] D. Grumiller and P. Romatschke, JHEP 0808, 027 

(2008). 

[27] J. L. Albaccte, Y. V. Kovchegov and A. Taliotis, 
JHEP 0807, 074 (2008), JHEP 0807, 100 (2008). 

[28] E. Avsar, E. lancu, L. McLerran and D. N. Triantafyl- 
lopoulos, JHEP 0911, 105 (2009). 



